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THE  MATHEMATICAL  BASIS  OF  EXPOSURE 
CONTROL  CALCULATIONS 


Introduction: 


The  current  exposure  control  system  Is  based  on  empirical  data  for 
standard  decay  which  is  shown  in  graphic  form  In  Attachment  A.  The  graphs 
indicate  a  linear  relationship* between  1og(R/Rj)  and  log  t: 

1. *  log^/Rj)  =  -1.2  log  t 

where  t  is  the  time  after  explosion  (in  hours)  of  the  radiation  reading 
Rt  (in  rads/hour),  and  Rj  is  the  radiation  reading  at  1  hour  after  the 

explosion.^ Note  that  this  is  standard  decay  with  a  decay  exponent  of  -1.2. 

» 

The  current  exposure  control  system  cannot  handle  non-standard  decay. 

/\ 

which  will  be  addressed  later. 

After  rearrangement,  equation  1  becomes: 

2.  Rt  -  Rj  t'1'2 

Integrating  with  respect  to  t  between  times  t  and  t.  (t  *  t.  )  gives  the 

a  D  a  D 

dose,  D,  received  during  that  time  interval: 

ftb  ,  o  ^b  ,  o  1 


D  =  /Rtdt 


t.  ft. 

“  -i  2  ^  .1  ? 

R.t  1  at  =  R.  t 
t.  ^t 


R1(-l/0.2)  -0.2t"1,2dt 


*  R1(-5)(t'°*2)|*b  =  -5Rj(tt 
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D  =  5R. (t 
1  * 


-0.2  *  -0.2, 


*An  index  to  the  variables  used  in  this  paper  (by  equation  number  and  page) 
is  included  as  Attachment  B. 
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Standard  Decay,  Current  Stay-Time: 


Stay-time  Is  based  upon  the  Maximum  Allowable  Dose,  MAD,  the  Total 
Accumulated  Dose,  TACD,  and  the  Allowable  Dose,  ALD,  for  a  trip  outside 


the  shelter: 


MAD  »  TACD  +  ALD,  or 

5.  ALD  *  MAD  -  TACD 
Using  equation  4: 

ALD  =  D  =  SRjU^0*2  -  tb"0,2),  or 

6.  ALD  =  5R1Ct"°'2  -  (t  +  CST)'0,23 

where  CST  is  the  Current  Stay-Time.  This  is  then  solved  for  CST: 
ALD/SRj  =  t-0*2  -  (t  +  CST)"0*2 


(t  +  CST)-0*2  =  t-0*2  -  ALD/5Rj 

t  +  CST  -  (t-0*2  -  ALD/SRj)-5 

7.  CST  =  {t-0*2  -  ALD/SRj)-5  -  t 

Solving  equation  2  for  R.  gives: 

1 

8.  Rj  =  Rt  t1'2 
Combining  equations  7  and  8  gives: 

9.  CST  =  (t-0*2  -  ALD/5Rt  t1*2)-5  -  t  ' 

=  t^"°*2^‘5)(l  -  ALD/5Rt  t)-5  -  t  jrj 

10.  CST  =  tC(l  -  ALD/5Rt  t)-5  -  1] 

Combining  equations  5  and  10  gives  the  standard  decay  CST  equation: 
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Standard  Decay,  Future  Stay-Time; 

Like  CST ,  Future  Stay-Time,  FST,  is  also  based  on  HAD,  TACD,  and  ALD, 
but  it  also  incorporates  the  dose  received  in  shelter  prior  to  the  antici¬ 
pated  exit,  D.  : 

in 


MAD  =  TACD  +  D.n  +  ALD,  or 


12. 


ALD  =  MAD  -  TACD  -  Dj 


Combining  equations  4  and  12  givesc 
13. 


ALD  =  MAD  -  TACD  -  SR^t-0*2  -  td"0,2)/PF 


where  t  is  the  current  time,  td  is  the  time  of  the  future  departure  (both 
in  hours  after  explosion),  and  PF  is  the  shelter's  protection  factor. 
Using  equation  4  again,  this  time  in  place  of  ALD,  gives: 

14. 


5R, 


Solving  for  FST: 


5RlCtd"°‘2  "  (td  +  FST)"°‘23  =  MD  -  TACD-^t'0,2-  td‘ 


td'0*2  -  (td  +  FST)'0-2  =  MAD5~  TACD  -  (t'°-2  -  td'0,2)/PF 
(tH  +  FST)'0,2  =  tH'0,2  -  ^AP-s-~- -TACD  +(t~0,2  -  t/1'°’2)/PF 


td  ♦  FST  = 


15. 


FST  = 


FST 


’  -0.2  MAD  -  TACD  ,  1  ,.-0.2  *  -0.2/ 

rd  '  5R.  PF  [Z  '  *d  } 

1  *-0.2  *  -0.21-5 
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*-0.2  *  -0.2 

*  -0.2  MAD  -  TACD  *  1  "  rd 

ld  '  5Rj  PF 

*  -0.2  .-0.2 

*  -0.2  .  TACD  -  HAD  .  "  rd 

*d  +  5Rj  +  - PF - 


-  t 
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Combining  equations  8  and  15  gives  the  standard  decay  FST  equation: 

,-5 


16. 


FST 


std 


*  -0.2  *-0.2 

*  -0.2  .  TACD  -  MAD  *  z  '  zd 
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Non-standard  Decay: 

One  major  flaw  In  the  current  exposure  control  system  is  its  inability 
to  allow  for  variable  or  non-standard  decay.  In  this  type  of  decay,  the 
"1.2"  exponent  of  equation  2  (and  subsequent  calculations)  varies  from  0.9 
to  2.0  as  a  result  of  weapons  variations,  overlapping  fallout  patterns 
(resulting  in  multiple  peaks  of  radiation),  weathering,  etc.  Realistically, 
this  must  be  taken  into  consideration.  (See  Effects  of  nuclear  Weapons , 
1977,  para  9.151.)  This  section  deals  with  the  derivation  of  CST  and  FST 
equations  for  variable  decay  exponents. 

Equation  2  becomes: 

17.  Rt  -  Rj  f" 


where  n  is  the  decay  exponent.  Therefore: 


18. 


0  =  /R tdt  * 


f*b 


R,t  dt  -  R 


'1 


hfry 


h. 


t~ndt 


(1  -  n)  t'n  dt 
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Non-Standard  Decay,  Current  Stay-Tiro; 

For  calculating  variable  or  non-standard  decay  CST: 

20.  ALO  *  MAD  -  TACD 

as  before;  then  using  equation  19: 

R, 


Solving  for  CST: 


'1 

=  [t1_n  -  (t  +  CST)1"0 


AL,D(n.  -_1).  -  _  (t  +  cST)1-n 

R1 

(t  +  CST)1'0  =  t1_n  -  ALD(n  -  1)/Rj 


t  +  CST  =  t1_n  -  ALD(n  -  1)/Rjj 

1 

CST  =  ft1-n  -  ALD(n  -  D/R.]1"0  -  t 


1/ ( 1-n ) 


Solving  equation  17  for  Rj  gives: 

23.  Rj  ■  Rt  tn 

Substituting  this  and  equation  20  into  equation  22  gives  the  variable 
(non-standard)  decay  CST  equation: 


CST  =  t1_n  + 
var 


*M. 
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Calculation  of  n  and  t: 

The  use  of  variable  decay  necessitates  a  method  of  evaluating  n  in 
terms  of  known  quantities.  Different  approaches  are  possible  depending 
on  whether  the  time  of  explosion,  tQ,  is  known. 

Known  time  of  explosion: 


| 

1 

i 

■  * 

i 

-I 


If  t  is  known,  data  from  two  radiation  readings  can  be  used  to 

« 

evaluate  n.  From  equation  17,  it  can  be  seen  that: 

R1  =  Ra  ‘a"  ‘  Rb  tb"  ’  s0 
30.  (yta>n  -  Ra/Rb 

.  n  log(tb/ta)  ■  log(Ra/Rb) 

;  )°9<w 
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Unknown  time  of  explosion: 

The  more  likely  case  would  be  that  tQ  is  unknown  or  that  several 
fallout  generating  blasts  occurred  at  times  different  enough  to  preclude 
the  selection  of  any  single  one  as  "the"  tQ.  In  this  case,  tQ  (or  an 
apparent  tQ)  must  be  calculated. 

Determination  of  n  in  this  case  requres  a  numerical  approach  and 
will  be  discussed  later;  for  now  we  will  assume  that  we  know  the  value  of 
The  derivation  of  t  is  based  upon  equation  30  and  the  following 
variable  notation: 

ta  =  time  elapsed  between  explosion  and  reading  a 

a 

tab=  time  elapsed  between  readings  a  and  b 

tfa  =  t  +  t  ^  =  time  elapsed  between  explosion  and  reading  b 

Inserting  the  last  equation  into  equation  30  gives: 

^b  _  la  +  fcab  .  /R  »l/n 
t  "  t_  '  iRa/Rb} 

a  a 


n. 
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The  Zulu  or  Local  time  for  the  time  of  explosion  can  then  be  determined 
from  ta  and  the  time  (Z  or  L)  of  reading  a. 

d 

f 

Calculation  of  n: 

To  find  n  in  terms  of  measurable  quantities  requires  a  third 
radiation  reading  at  time  tc.  Using  equation  32  and  the  variable  tbc 
to  represent  the  time  between  readings  b  and  c: 


>v 


Simplification  and  multiplication  by  1  gives: 

he  __  (ync)1/n  -  .  (yy1/n 

*ab  <Ra/Rb)1/n  -  1  (VRa)1/n 

tbc  (yRc)1/n  - 1 
*ab  1  -  (Rb/Ra)1/n 

Choosing  the  time  intervals  between  the  radiation  readings  so  that  tbc = tgb 
gives  an  .intuitively  interesting  result: 

,  ‘be  (VRc>1/n-1 
tab  1  -  (Rb/Ra)1/n 

1  “  ^Rb/Ra)1/n  =  <Rb/Rc)1/n  _  1 

(Rb/Rc)1/n  +  (Rb/Ra)1/n  =  2 
Vn  (R  -1/n  +  R  -1/nj  .  2 

DC  a 


34. 


R  "1/n 
Rb 


Note  that  equation  34  resembles  an  averaging  formula.  Also  note  that  this 
equation  does  not  easily  transform  into  an  explicit  equation  for  n; 
therefore,  we  will  resort  to  an  approximation  method  to  find  n.  Continuing 
from  equation  33: 

‘bet1  -  <VRa>1/n)  ■  ‘ab«VRc>1/n  - 

‘be  *  ‘be<VRa>l/"  -  ‘ab<VRc>1/n  '  ‘ab 
‘ab  +  ‘be  *  ‘ab<Rb/Re>1/n  ♦  ‘be'W1'" 


T  =  tab(Rb/fV  +  WW 

35.  f(n)  =  1  =  ^  (Rb/Rc)1/n  +  ~T~  (Rb/Ra)1/n 

R  1/n 

36.  f(n)  =  1  =  -V  ^ab^RC1/n  +  * be /Ra1/n> 

where  T  =  +  tbc.  Equation  36  is  an  implicit  equation  for  n  which  can 

be  used  to  develop  a  numerical  approximation  algorithm  by  which  n  can  be 
calculated  to  any  desired  level  of  precision  on  a  programmable  calculator. 
Before  developing  the  algorithm,  the  theory  behind  the  approximation  method 
will  be  discussed. 

Newton's  Method  of  Approximation: 


For  any  function  y  =  f(x)  at  point  (w^,  f (w^ ) ) ,  the  slope,  m,  of  the 
tangent  line  is  given  by  the  equation 


m  =  f'(w,)  =  — — 


f(x)  -  f(w.) 


wV.V.V.  V*  VLV_\  S.  _V V*  -&.7JV  N  -'•  .v  .  *  _  *  :  •  .  ' 


where  f(x)  is  known  and  is  taken  as  the  y-value  on  the  tangent  line  at  w2, 
and  Wj  is  the  first  "best-guess"  for  an  approximate  value  for  x.  Solving 
for  w2  gives: 

f(x)  -  f(Wj) 

W2  =  W1  +  ■  f'(*p - 

Note  in  the  figure  on  the  preceding  page  that  w2  more  closely  approximates 
x  than  does  Wj,  and  that  repeated  application  of  the  procedure  ing  w2 
in  place  of  Wj  results  in  w^,  an  even  better  approximation  for  Therefore, 
the  general  approximation  equation  is: 


x  ~  w  +  f(x).  ~ 

x  w  f '  w) 


where  w  is  the  current  best  approximation  for  x.  The  procedure  can  be 
iterated  until  an  adequately  precise  approximation  is  reached. 
Approximation  of  n: 

Equation  38,  then,  is  the  basis  for  estimating  n: 

38.  n  =  w  +  Vll|wl 

where  w  is  an  approximation  for  n,  f(x)  =  f(n)  =  1  (from  equation  36),  and 
f(w)  is  calculated  using  w  for  n  in  equation  36.  Next,  f'(w),  the  first 
derivative  at  n  =  w  needs  to  be  calculated.  Starting  from  equation  35: 

f(n)  =  V  (Rb/Rc,n"1  +  ¥  'W"'1 

f '  (n)  =  ^  (Rb/Rc)n  1("n"2)  ln  (R b /RC)  + 

¥  (vRa)n"1(-n_2nn  (Rb/Ra) 


f '  (n)  = 


-T7T  In  (Ru/Rj  +  — ~ttt  In  (Rk/Rj 


To  find  the  best  approximation  for  n,  the  equation  is  used  iteratively, 
taking  the  resultant  n  from  each  step  for  the  value  of  w  in  the  next  step  until 
the  desired  level  of  precision  is  reached.  A  reasonable  limit  in  this  case  is 
43.  jn  -  w|  <  0.05 


This  method  of  approximation  rapidly  converges  on  n  so  that  the  limit 
should  be  reached  within  a  few  iterations.  When  n  is  found,  it  can  then 
be  used  with  equation  32  to  find  the  time  in  hours  since  the  explosion. 


This  can  then  be  used  with  n  and  the  current  radiation  reading  in  equations 
24  and  29  to  find  the  Current  and  Future  Stay-Times.  (An  approximation 
procedure  similar  to  the  one  described  here  could  be  used  with  equation  29 
to  calculate  the  time  of  departure  from  shelter,  td,  at  which  FST  will  equal 
the  amount  of  time  projected  to  be  spent  outside.  In  other  words,  an  opum*l 
departure  time  could  be  calculated.) 

Conclusion: 

It  is  important  to  note  that  using  this  analysis,  exposure  control 
calculations  require  knowledge  of  only  three  consecutive  (and  decreasing) 
radiation  readings  and  the  two  time  intervals  between  the  readings.  For 
this  set  of  equations,  there  are  no  restrictions  on  the  time  intervals. 

Most  importantly,  the  analysis  presented  here  does  not  assume  standard  decay, 
nor  does  it  require  knowledge  of  the  time  of  the  nuclear  detonation(s) . 


Reference:  The  Effects  of  Nuclear  Weapons,  Samuel  Glasstone  and  Philip  J. 
Dolan,  eds.,  1977.  (Dept,  of  the  Army  Pamphlet  50-3). 


Revised  1  March  1983 
S.  Mark  Clardy 
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VARIABLE  INDEX 
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1 

Variable 

Definition 

First  Appearance 
Eq 1 n  #  Page 

A 

tbc/Ra^W  (substitution  variable) 

41 

12 

ss 

ALD 

Allowable  dose  (rads) 

5 

2 

a 

B 

T/Rb^W  (substitution  variable) 

41 

12 

•> 

C 

tab/Rc1/W  (substitution  variable) 

41 

12 

CST 

Current  stay- time  (hrs) 

6 

2 

a 

D 

Dose  (rads) 

3 

1 

S 

l' 

Din 

Dose  received  in  shelter  before  exit  (rads) 

12 

3 

FST 

Future  stay-time  (hrs) 

14 

3 

MAD 

Maximum  allowable  dose  (rads) 

5 

2 

n 

Decay  exponent 

17 

4 

v- 

PF 

•Protection  factor  of  shelter  (=  j^side6  p) 

13 

3 

R 

Radiation  reading,  current  (rads/hr) 

11 

2 

* 

R1 

"  "  at  1  hr  after  explosion 

1 

1 

a 

§ 

"a 

"  "  "  time  ta 

a 

7 

Rb 

II  II  II  II  A. 

zb 

7 

Rc 

M  II  It  II  a. 

C 

9 

r 

Rt 

"  "  "  "  t  (indefinite  time) 

1 

1 

T 

tab  +  *bc  (substitution  variable) 

35 

-  4 

t 

Time  after  explosion  (hrs) 

1 

1 

4a 

Time  elapsed  between  explosion  and  R,  (hrs) 

a 

32 

8 

■> 

4b 

n  n  ii  it  n  n 

Rb 

- 

8 

•i 

4c 

ii  n  ii  ii  ii  p 

- 

8 

4d 

"  "  "  "  "  departure 

from  shelter 

13 

3 

*ab 

"  "  "  readings  Rg  and  Rb 

32 

8 

4 

tbc 

"  "  "  "  Rh  and  R_ 

b  c 

9 

r- 

w 

Approximation  for  n 

38 

11 

14 


